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We describe a plausible-speculative form of quantum computation which exploits particle 
(fermionic, bosonic) statistics, under a generalized, counterfactual interpretation thereof. In the 
idealized situation of an isolated system, it seems that this form of computation yields to NP- 
complete=P. 



I. INTRODUCTION 



Particle (fermionic, bosonic) statistics has never been applied until now in the usual (most fruitful) approach to 
I quantum computation. We mean sequential computation, which can be characterized by the fact that a reversible 
Boolean network appears in the time-diagram of the computation process^ 1-6 '. 

It has been applied instead in the approach expounded in ref. [7-12], which differs from the usual one: 

> ' 

• a different form of computation is used; its classical counterpart is simulated annealing or ground state compu- 
tation; time is orthogonal to the layout of a reversible Boolean network whose input and output qubits coexist 
as simultaneous eigenvalues of a set of compatible observables. Therefore part of the network input and part of 



OO 
O 

' the output can be both constrained. Checking whether this network is satisfiable is a version of the NP-complete 

00 ! SAT problem. 



This framework seems suitable to applying particle statistics: the network gates and wires establish Boolean rela- 
tions between the simultaneous eigenvalues of the respective input/output qubits; a particle statistics symmetry 
is also a constraint applying to simultaneous eigenvalues. 



• However, in the conventional interpretation, a particle statistics symmetry is an initial condition which is 



conserved as a constant of motion by a unitary evolution. Therefore it does nothing along the course of this 
type of evolution. In the generalized interpretation propounded in ref. [8-12], a particle statistics symmetry 
is the result of a continuous projection on the Hilbert subspace satisfying the symmetry. We can thus apply 
a Hamiltonian G r to a part r of the system, provided that particle indistinguishability is preserved. This and 
the continuous projection, seen as a form of "interaction" between r and the remaining part of the system s, 
define a system level Hamiltonian G rs . G rs induces a unitary evolution of the state of the whole system which 
does not violate the symmetry. In a way, this evolution is driven by G r and shaped by particle statistics seen 
as continuous projection. 

In [12] we have shown how to map the logical constraint established by a network element (N.E., namely gate or 
wire), on a constraint induced by particle statistics. This constraint is seen as continuous projection on the Hilbert 
subspace satisfying the N.E. logical constraint. The network is prepared in an initial state satisfying the (partial) 
input constraints and all N.E. constraints, whereas the output constraint (it is sufficient to constrain only one output 
qubit) is temporarily removed. This amounts to solve a problem polynomial in network size. Then we operate on 
the (to be constrained) output qubit in order to bring it in match with its constraint. This operation, under the 
continuous projection, induces a unitary transformation leading the network state to satisfying all the constraints 
(with very high probability), if the network is satisfiable. This requires a time independent of network size and would 
yield to NP-complete=P, see [12]. 

However, in [12] we could not exclude the possibility that the operation performed on the output qubit pushed the 
network into an error state, with some probability; the rate of growth of this probability with network size was not 
known. An exponential growth would have completely vanificd the former result. 

In the current work, we show that this probability can in principle be zero, which brings to NP-complete=P in an 
idealized framework. 
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To show this, we will need to recover part of [12], namely the counterfactual interpretation of particle statistics and 
its application to an elementary network element (a NOT gate). 



II. A COUNTERFACTUAL, GENERALIZED INTERPRETATION OF PARTICLE STATISTICS 

The first part of this Section summarizes a corresponding part of [12]. The second part introduces a further 
development. 

Let us consider the unitary evolution of a triplet state induced, for assumption, by particle statistics: 

|* (t)) = cos 2 (i? + uot) |0) x |0} 2 + sin (tf + uot) cos (tf + uot) (|0) 1 |1) 2 + |0) 2 ) (1) 
+ sin 2 (# + u>t)\l) 1 \l) 2 

1 and 2 label two identical two-state particles with Hilbert spaces respectively Hi = span{\0) 1 , |1) 1 } and 7Y 2 = 
span{\0} 2 , |1) 2 } (which makes this an idealized situation). 

Evolution (1) is obtained by applying an identical rotation uot of the states of the two particles to the initial state 

|* (0)), namely: Q 1 (uot) Q 2 (uot) |tt (0)) = |* (t)) , where Q j (uot) = (_ C °^ f™Jt ) ' with j = 2 " 

In the usual interpretation, the propagator Qi (uot) Q 2 (ut), symmetrical for the particle permutation Pi 2 , commutes 
with S12 = \ (1 + P12)', therefore the triplet state symmetry is an initial condition conserved as a constant of motion. 
The counterfactual interpretation given in [12] can be summarized as follows. We consider the possibility that 
(t)) at any time r = t — dt goes out of symmetry. This is interpreted by saying that |* (t — dt)) is not constrained 
by the symmetry: it is therefore a free vector of the Hilbert space H12 = Hi (g> H2, namely: 

1* (t dt)} = 4~ dt) io)! io) 2 + \ (4i dt) + ^- dt) ) mi ii) 2 + io) 2 ) + (2) 

\ (<~ dt) - at dt) ) (\0)i |i) 2 - 11), |0) 2 ) + «<«"*> \i)i |i) 2 ; 



the a£* dt ^ are free complex variables independent of each other up to je{o 1} 
other words, if ti 7^ t 2 , \^ (ti)) and |\& (t 2 )} are two independent free normalized vectors of 7ii 2 . 

Of course we must assume the possibility that a^ dt " > — a^*' 7^ 0. If this is the case, namely if |* (t — dt)) is out 
of symmetry, this state should be "immediately" projected on the symmetric subspace 

W# = span jlOK |0) 2 , (10), |1) 2 + |l) x |0) 2 ) ,|1) 1 |1> 2 

The result of this projection can be obtained by submitting |^ (t)), another free normalized vector of H12 indepen- 
dent of |^ (t — dt)), to the following equations: 
for all t: 

i) Si 2 |*(t)) = |*(t)>, 

ii) Max |(*(i) | &(t — dt))\ , which naturally means that the distance between the vector before projection 
I'J' (t — dt)) and the vector after projection |^ (t)) is minimum, 

where |* (t)) is a free normalized vector of Hi2- 

These equations also take into account the fact that projection goes on repeatedly in a continuous fashion. 

In this way, the triplet state symmetry is seen as the result of a watchdog effect that continuously projects the 
system state on the symmetrical subspace • Of course the result is that |* (t)) can never go out of symmetry. 
However, this counterfactual reasoning will have consequences (the possibility that counterfactuals yields to effective 
consequences in the quantum framework has been highlighted by R. Penrose^ 13 !). 

A first consequence is that we do not need to assume that symmetry is an initial condition: no state satisfying (i) 
for all t can ever be out of symmetry. 

Also the notion that symmetry is a constant of motion can be given up. To show this, it is useful to see the 
continuous projection (i) and (ii) as a special, continuous form of state vector reduction. As a matter of fact, the 
cancelation of the amplitude aoi _d ^ — ^io ^ an d the renormalization of the other amplitudes in state (2) is a partial 



(t-dt) 2 
ot - ■ 
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1, for any time t. In 
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state vector reduction on the subspace 7i\ 2 (it is partial since this subspace has dimension higher than one). Of 
course there is no dynamics in state vector reduction, only interference - in the above form of cancelation of one 
amplitude and renormalization of the others. 

We should note the peculiarity that such a reduction always occurs on the symmetrical subspace, never on the 
orthogonal subspace. This is paradoxical in counterfactual reasoning, it would mean that the result of reduction is 
not random but is affected by a condition placed in the immediate future, namely that the result does not violate the 
symmetry. As a matter of fact, the current counterfactual interpretation can be justified in a two-way (advanced and 
retarded) propagation model ^ . 

From another standpoint, there is no paradox at all since, actually, 1^ (t)) never goes out of the subspace Ti.\ s 2 . 
However, the consequences of this counterfactual reasoning will diverge from the conventional way of applying quantum 
mechanics. 

We shall now reconstruct evolution (1) by resorting to conditions (i) and (ii). A first observation is that these 
conditions affect the overall state of the two particles, therefore the states of the individual particles may no more be 
defined. We must use the particle density matrices 

Pi (t) = Tra_i[|* (<)><* (t)|], 

where Tr^-i means partial trace over 3 — i, and i = 1, 2 is the particle label (if we use the method of random phases' 1 ^ , 
| "J (t)) does not need to be a pure state - anyhow it will turn out to be that). 

A second observation is that the coherence elements of each density matrix - as entanglement - can also be affected 
by the watchdog effect (in fact they will be determined by it). We only know that the diagonal of each density matrix 
must show an u>t rotation [to be consistent with eq. (1)]. In conclusion we must add to conditions (i) and (ii) the 
further double condition (iii): 

diag Pl (t) = diag {Tr 2 [|* (t)) (¥ (t)\}} = cos 2 (tf + cot) |0) x <0| x + sin 2 (0 + cot) (l^ . 



diagp 2 (t) = diag{T ri [|* (t)> (tf {t)\}} 



{d + cot) |0) 2 (0| 2 + sin 2 (0 + cot) |1) 2 (1| 2 



It is readily seen that the simultaneous application of conditions (i), (ii) and (iii) yields evolution (1). It can be 
said that evolution (1) is driven by conditions (iii) and shaped by conditions (i) and (ii)Q 

We should now note a fact which is essential to the current work. As readily seen, by removing either one of the two 
conditions (iii), evolution (1) is still obtained - the two conditions are redundant with respect to one another. 

Therefore it is perfectly legitimate to say that the rotation of the state of only one particle, either one in an 
indistinguishable way, drags an identical rotation of the state of the other particle. In this idealized picture, particle 
statistics can be seen as an interaction free constraint, namely as a non-dynamic constraint operating by way of 
destructive interference and renormalization. 

We will show how to apply (speculatively) the driving condition (iii). This is a further development with respect 
to [12]. Only one particle j, either one in an indistinguishable way, should be submitted to the Hamiltonian 



Gj = LOO 



(3) 



where a y j is the Pauli matrix in 7ij 



i 
-i 



(4) 



and la—j is the identity in Tia—j. The overall Hamiltonian G12, operating in H.12, is obtained through symmetrization 
of Gj-. G12 = \cj (o-yila + lic^)- This yields 



G12 = -LJ 



( i i 0" 
-i i 
-i i 
V -i -i 0, 



(5) 



1 An example of a unitary evolution shaped by a continuous form of state vector reduction, is the evolution of the polarization 
of a photon going through an infinite series of polarizing filters, each rotated by an infinitesimal constant amount with respect 
to the former one. In a way, we go back to the root of quantum computation (computation reversibility' 15,16 ') and take an 
alternative branch, by exploring a strictly quantum form of reversible computation. 
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This is in fact the generator of Q\ (tot) Q2 (u>t), as readily checked. 

The above symmetrization is interpreted as follows. If we applied the Hamiltonian Gj to one particle assumed to 
be independent of the other, we would have obtained a rotation of the state of that particle by tot. By considering the 

[si 

continuous projection on Ti\ 2 as a form of interdependence between the two particles, the former Hamiltonian becomes 
G\i at overall system level (in a way, continuous projection is considered as a sort of interaction Hamiltonian between 
the two particles). G\i generates the operator Q\ (cot) Q2 (u>t) which rotates the states of both particles of the same 
amount. Therefore we can say that the rotation of the state of only one particle (either one in an indistinguishable 
way), drags an identical rotation of the state of the other. 

This, in the current context, is a tautological interpretation of particle statistics. However, applied to a different 
context (this will be two indistinguishable particles hosted by two distinguishable lattice sites), such an interpretation 
will yield two results diverging from the conventional way of applying quantum mechanics. 



III. AN ELEMENTARY GATE AS A PROJECTOR 



We consider a couple of (coexisting) qubits r and s which make up the input and the output of a NOT gate. Let 

H rs = span{\0) r |0) s , |0> r |l) s , |l) r |0> 8 , |l) r |1> J , 
be the Hilbert space of the two qubits, 

Tig =apon{|0) r |l) 8 ,|l) r |0)J 

be the constrained subspace, spanned by those Ti rs basis vectors which satisfy the NOT gate, and A rs be the projector 
from H rs on Tig ■ 

We shall apply the mathematical model of Section II to represent an evolution of the state of the NOT gate \^ (t)} 
(this Section is purely mathematical, a plausible physical model will be given in Section IV). (t)) should be a free 
normalized vector of Ti rs , 

!*(*)>= E ^V)r\3)s, With ]T |ag )2 
i,j£{0,l} »,j€{0,l} 

subject to continuous projection on Tiers'- 



for all t: 



i) A rs \*(t)) = |*(i)>, 

ii) Max (i) \*(t-dt))\. 

Let I W (0)) = cosfl |0) r |l) s + sin$ |l) r |0) s be the gate initial state. We assume of acting on qubit r with the 
Hamiltonian 



G r = LU(Ty r , where a yr = lu 




We have two ways of deriving the evolution induced by G r . 

a) If acting on the independent qubit r, G r would rotate its state by tot. As we have seen in Section II, the diagonal 
of qubit r density matrix should not be affected by the continuous projection on T6?} '. Thus: 

hi) diag Pr (t) = diag {Tr s (t)) (* [t)\]} = cos 2 (d + cut) |0) r (0| r + sin 2 (d + cut) |l) r (l| r , 

Conditions (i), (ii) and (hi) define the unitary evolution 

I* (t)> - cos (<? + cut) |0) r |l) s + sin (0 + cjt) |l) r |0) s , (6) 

as readily checked. Condition (hi) drives and conditions (i) and (ii) shape this evolution. We can see that the rotation 
of qubit r [the driving condition G r , or (hi)] induces an identical rotation of qubit s: 
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Ps (t)=T rr [|*(t)> <* 



= sin 



( l ? + ^)|0) s (0| s +cos 2 (0 + W t)|l>,<l| 



of course and 1 are interchanged. 

b) A second way of deriving evolution (6) consists in computing the Hamiltonian G rs acting on the overall state 
of the two qubits. G rs is originated by G r and the continuous projection (seen as a form of interaction or better 
interdependence between the two qubits). Given that this latter introduces the following bijective correspondence 
between H r and H rs - 



\0) r " |0) r |l) s , |l> r ~|l> r |0) a) 



G r is: 



/° 














i 





—i 





Vo 









naturally G r populates the central part of G rs . 
As readily checked, G rs generates the operator: 



Qrs (vt) 



with|0) r |l) s = 




One can see that Q rs (tot) |\& (0)) = (t)) given by equation (6). We have thus obtained the same result as before 
- when the driving condition was the evolution of diagp r (t). 

In view of what will follow, it is important to note that, if all 7ifs basis vectors occur with amplitudes different 
from zero in the initial state (t = 0), namely if $ ^ 0, condition (i) is redundant with respect to condition (ii). In 

this case, condition (ii) alone implies a^o = = which already satisfies condition (i). 

On the contrary, condition (i) is not redundant if i3 = or ^. For example, if $ = 0, i.e. (0)) = |0) r |l) s , 

2 



condition (ii) 



constraint 



,(*) 
no 



implies a 

2 



and 



(*) 
% 



cos 2 (i9 + wt) , as needed, while a^J and are only subject to the 



(*) 



forbidden state a'^ |l) r |l) s . 

In conclusion, we have ascertained a peculiar fact. Our "operation on a part" [this is just the mathematical 
condition (hi), or G rs in equivalent terms, for the time being], blind to its effect on the whole, performed together 
with continuous A rs projection, generates a unitary transformation which is, so to speak, wise to the whole state, to 
how it should be transformed without ever violating A rs (i.e. the NOT gate). Of course A rs ends up commuting with 
the resulting overall unitary propagator, but because this is itself shaped by A rs . 



sin 2 (# + cot) . Thus, disregarding condition (i) would allow for the existence of the 



IV. EXPLOITING PARTICLE STATISTICS 



A rs projection can be shown to be an epiphenomenon of particle (fermionic or bosonic) statistics "turned on" in a 
special physical situation. In the following, we will adopt fermionic statistics. 

In order to implement the NOT gate, we consider two identical fermionic particles 1 and 2. Just for the sake of 
visualization, we can think that each particle has spin 1/2 and can occupy either one of two distinguishable lattice 
sites r and s. The generic basis vector of this system has the form |xi)i 1X2)2 l^i)i | A 2 ) 2 7 where %i = 0, 1 means 
that the spin of particle i is down, up and \ = r,s means that the site occupied by particle i is r, s. For example, 
|0) 1 |1) 2 |0i l s )2 rea ds: particle 1 spin = 0, particle 2 spin = 1, particle 1 site = r, particle 2 site = s. 

16 combinations like this make up the basis of the Hilbert space Hi2- However, there are only six antisymmetrical 
combinations (not violating fermion statistics) which make up the basis of the antisymmetrical subspace ■ 

These basis vectors are represented in second quantization and, when there is exactly one particle per site, in qubit 
notation (x and A stand respectively for the qubit eigenvalue and label), |0) is the vacuum vector: 

I a) = a\ r a\ r |0) , 

|6) - al a\ s |0> ; 
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|c) = al al |0) - |0) r |0) s , 
\d)=a\ r aljO) = |l) r |l) s , 

l e ) = 72 ( a l 4s + 4r 4s) 10) = ^ (|0) r |l) s + |l) r |0) s ) . 

If) = 7= 2 (4r 4s - 4r 4s) |0) = ^ (|0) r |l) s - |l) r |0) J . 

a\^ x creates a particle of spin \ m site A; creation/annihilation operators are subject to: |a^,aj| = {a,i,cij} = 

0, |a- , a^j = 5i t j. Under the condition that there is exactly one particle per site, they generate a qubit algebra. 
Now we introduce the Hamiltonian 

H rs = —{E a a 0r a\ r a 0r a lr + E b aj s a\ s a 0s a ls + E c a) 0r a 0s aor a o« + E d a\ r a ls a\ r ai s ), 
with E a , Eb, E c , Ed> E discretely above 0. This leaves us with two degenerate ground eigenstatcs: 

|e) = (|0) r |l) s + |l) r |0) s ) and |/) = ± (|0) r |l) s - |l) r |0> J . 

The generic ground state is thus: 

l*> = <* I0) r |l) s +P\l) r \0). , with M 2 + |/3| 2 = 1. (7) 

Of course |\&) satisfies A rs |*) = \^>), and belongs to TL r c } (Section III), a subspace of Ti[^ . 

Let A\2 I 1 ?) = \ (1 — Pyi) be the usual antisymmetrization projector. Due to the above anticommutation relations: 

M2 |0> r |l) s = |0) r |l) s and A 12 |l) r |0) a - |l) r |0) a , moreover 

^12 |0) r |0} s = |0) r |0) s and A 12 |l) r |l) s = |l) r |l) s , 

without forgetting that |0) r |0) s = |c) and |l) r |l) g = \d) are excited states. 

The NOT gate can be implemented by suitably operating on the ground state (7). We assume this to be initially: 

|*(0))=cosi?|0) r |l) s +sint?|l} r |0) s . (8) 

The evolution induced by the Hamiltonian G r is obtained by adopting the projection interpretation of particle 
statistics. This means that \*5? (t)} is continuously projected on the antisymmetric subspace while diag p r (t) 
evolves as if G r were applied to the independent qubit r: 
for all t: 

i) ^ 12 |*(t)) = |*(t)), 

ii) Max (i)| *(t-dt))| , 

iii) diag p r (t) = diag{Tr s [(|* (t)) (* (t)|)]} = 
cos 2 (0 + wi) |0) r (0| r + sin 2 (i? + wt) |l) r (l| r , 

iv) (£„(*)) = (*(t)|IT„|*(t)>=0 ) 

where |* (i)} is a free normalized vector of Hi2- 

The solution of the above equations is the desired evolution (6), repeated here for convenience: 

|* (t)) = cos (i? + ut) |0} r |l) s + sin (0 + wt) |l) r |0) s . 

Conditions (i), (ii) and (iii) mean that the link state undergoes a transformation [driven by (iii) or G r ] under continuous 
state vector reduction on the antisymmetric subspace Ti.[^ ■ If i? ^ 0, ^, namely if the preparation (4) comprises all 

the basis vectors of H r °J , condition (ii) alone keeps the link evolution inside H r c ) (Section III) . The link state remains 
ground and consequently the link expected energy (£ rs (t)) is always zero. By excluding # = 0, |, condition (iv) is a 
consequence of the former conditions. 

Mathematically, conditions (i) and (iv) give the constraint A rs |*E r (t)) = \ty (t)). In conclusion the above conditions 
(i) through (iv) (which imply interpreting fcrmionic antisymmetry A12 as continuous projection on the antisymmetric 
subspace) are equivalent to condition (i) through (iii) of Section III. This gives in fact the evolution (6). 

Since [G rs , H rs ] = 0, as readily checked, the application of G r to qubit r docs not disturb in principle the ground 
state of H rs , in constrast with what hypothesized in ref. [12]. This means NP-complete=p[ 12 l under the current form 
of computation. Of course we should keep in mind that we are in the same idealized and speculative context of the 
former work. 
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V. CONCLUSIONS 



By using the generalized interpretation of particle statistics (viewed as continuous projection on a constrained 
subspace), we have obtained an evolution of the NOT gate by acting only on qubit r (an input or an output qubit of 
the gate, indifferently so given that the two coexist) and without violating the NOT gate logical constraint. 

This would be impossible under the conventional interpretation of particle statistics. In this context, since all TL rs 
basis vectors are already antisymmetrical, particle statistics would do nothing. Any operation performed only on 
qubit r in the initial entangled state (8), would necessarily originate the terms |0) r |0) s and |l) r |l) s , thus violating 
the NOT gate logical constraint. 

If the model Hamiltonians used in this work could be substituted by more concrete Hamiltonians - say implementable 
in a laboratory - in principle such a divergence between the two interpretations of particle statistics could be verified. 

This research has been partly developed during the Elsag Bailey-ISI Workshop on Quantum Computation (Turin, 
1997). Thanks are due to A. Ekert, D. Finkelstein, S. Lloyd and V. Vedral for useful suggestions. 
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